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ABSTRACT 

Although organized transportation systems are very old indeed, 
active research into their operations is relatively new. Some of the 
problems encountered in transportation research are discussed along 
with methods for their solution, The Hitchcock distribution problem 
is discussed and a sample problem solved using three different methods. 
The usefulness of graph theory in Studying transportation problems is 
pointed out. Examples are then given of fits use in a maximal network 
flow problem and in the problem of minimizing equipment requirements 
to meet a fixed schedule. Finally, some miscellaneous techniques are 
discussed, 

The writer wishes to express his appreciation for the guidance and 
encouragement given him by Professor Thomas E. Oberbeck of the U. S. 
Naval Postgraduate School, for the helpful suggestions of the second 
reader, Professor H. C. Ayers, and for the careful clerical assistance 


provided by Mrs. Norma A. Stevens. 
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INTRODUCTION 

Im the past decade research into tramsoorcation problems has been 
greatly imtensified. Much of this research has involwed the application 
of the principles and techniques from 4 wide variety of mathematical 
disciplines. In this thesis we shall summarize some of the existing 
techniques and give examples of their applicability. 

Most of the research we have found reported has as its object the 
optimization of some feature of a transportation activity. Some examples 
are: minimizing costs, relieving congestion, minimizing equipment re- 
quirements, establishing economical maintenance schedules, minimizing 
loading and unloading times. Wsually the methods reported are directed 
at solving a problem of a particular transportation activity such as a 
railroad or an airline. But much of the research has been sufficiently 
general to be widely applicable. The tendency to zeneralize models is 
prevalent. For exampie, a researcher may seek to minimize the number 
of "units" required to maintain a fixed schedule, "Wnits'’ may then be 
interpreted by a railroad as boxcars or trains; by am airline, as air- 
planes; and by a shipping company, 45 tamkers or freighters, 

Transportation systems of any appreciable size gre extremely 
complex entities and may defy analysis im many respects, Researchers 
turn to mathematical models of the systems they wish to analyze and 
build into the model the degree of detail or generality necessary to 
give significant results. This practice has resulted im the application 
of a great many mathematical techniques to problems in transportation. 
The techniques and disciplines which have been used include: linear pro- 
gramming, transhipment and assignment techmiques, non-linear and dynamic 
programming, queuing, or waiting line theory, graph theory, and the Monte 


Carlo method, 1 





In our research for this paper we nave found 4a great deal of mater- 
ial on transportation type preblems in the following publications: The 
Journal of the Operations Research Society of America, Management Science, 
Naval Research Logistics Quarterly, and others listed in bibliographies 
at the end of each charter. Groups we have found to be particularly 
active in transportation research include’ The Rand Corporatiizn, Santa 
Monica, California; the George Washington University Logistics Research 
Project; The Cargo-Handling Research Project at the iniversity oe Calif- 
ornia at Los Angeles; The Management Sciences Research Group, Purdue 
University; to name only a few. 

Professor F, Harary of the University of Michigan has done a consider- 
able amount of research in graph theory, which is a directly supporting 
discipline, He is co-author of a new bock om this subject which is due 
to appear in 1960. This book will be the first full treatment of this 
subject in English. 

As an indication of the direction in which current research is tend- 
ing, we point to Rand report R-351 which gives, im abstract form, the pro- 
ceedings of "The Rand Symposium on Mathematical Programming’’. In his 
address to the opening session of this symposium, G. B. Dantzig cited the 
five following theoretical areas as those in which future developments 
are expected; 

(1) Special Structures 
(2) Discrete Programming 
(3) Network Theory 
(4) Non-linear Programming 
(3) Uncertainty 
It is anticipated that transportation theory will be one of 4a great 


many béneficiaries of developments in these areas. 











Tn arranging various types of transportation problems irto groups 
for the ensuing discussion wea were, of necessity, a bit artitrary. We 
have presented in detaii only a few of the more widely used cechniques 
and have made reference to others, We have tried to present a reason- 
ably complete bibliography at the end of each chapter. 

In Chapter I we have stated and discussed some methods for the 
solution of the Hitchcock distribution problem, both with and without 
capacity limitations. These are essentiully linear programming pro- 
blems. Im Chapter II we have indicated the applicability of Linear graph 
theory to transportation problems and general network problems. In 
Chapters III and IV we have we have discussed in some detail two problems 
that illustrate the applicability of graph theory to transportation 
problems. In Chapter V we have cited several examples sf other techniques 
and the problems to which they have been applied, Specifically, wa have 
discussed queuing theory, simulation, and straightforward engineering 
analysis including the difficult and time consuming task ot cate gather- 
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CHAPTER I 
THE TRANSPORTATION PROBLEM IN LINEAR PROGRAMMING 
1. The Hitchcock Distribution Problem. 

The problem to be discussed im this chapter was originally stated 
in 194i by F. L. Hitchcock Gel: and was further discussed and elaborated 
by T. C. Koopmans C2]. The problem may be stated as follows: Certain 
amounts, a, (i = 1,...,m), of a homogeneoug product are available at 


each of m origins, and certain amounts Be (j = 1,.-.00,m) are required 


at each of n destimations. The cost, ,» of shipping a unit amount 


Cis 
th ‘ th a : , 
from the i origin to the } destination ts known for all i and j. 


We will further stipulate 24, = Zb,. The problem 1s to select the 


amounts to be shipped form each origin to each destination so that 


*i4 


the total shipping cost is 4 minimum, The amounts aE must be non- 


i 


Algebraically, we state the problem in the tollowing form: 


negative, x fe Ge 


Given: (1) Two sets of cunstants ay df, aerate at 
and dD, j = pyrarate 0 


such that 


Sa. “Ze, and 


L 


: __ r aa 
(2) A matrix of constants C= « cy J t= Ly oso gt 


Find: A set of non-negative unknowns 


x ees OS satisfying 


ij 7 
a Ss x, = b-?> = 1i,4te,% 
@ Sx, = yg f=, 
(b) : > ae = a. — ites so, 0 
— ij i 
) 
1, Numbers in brac ets refer to references cited in the Bibliography 





eae 
(Cc) pas c., «,. = Minimum 
. ‘ ’ Lj tj 


~ ee 


Linear programming problems which may be stated in the foregoing 
form have come to be calied "the transportation problem” even though 
they do not all deal with transportation, Examples of such problems are 
the personnel assignment problem, the contract - award problem, [3 | 
{4} the Traveling-Salesman problem fio] , and the Caterer problem (1) 0 

G. B. Dantzig [5] was first to formulate this problem ee a linear - 
programming problem, and he gave an algorithm for its solution based on 
his Simplex method lé | » The transportation algorithm is considerably 
easier to apply than the general simplex method in that it does not re- 
quire the inversion of matrices. The algorithm is discussed in more 
detail and a sample problem worked out in Appendix A. 

2, The Assignment Technique. 

As was noted above, a special case of this "transportation pro- 
blem"’ is the n x n assignment problem. In this case all a, = on aaa 
and all sa are either one or zero. The problem may be stated as 
follows: determine the optimum assignment of nh persons to mn jobs, 
given a matrix of constants r, J which specify, in appropriate units, 
the worth or performance rating of the igh individual in the ae job. 
The most straightforward approach is to consider each of the ne: possible 
arrangements, But, for only moderately large n, this becomes a hope- 
less task, What is required is a procedure for finding an optimum assign- 
ment (there may be more than one) with a reasonable amount of effort. 


Such a procedure was provided by H. W. Kuhn [3] in his "Hungarian Method". 


Kuhn's algorithm has been modified and extended to the mk n 





transportation probiem by L. R. Ford. Jr. and D. R. Fulkerson [7] : 

and by J. Munkres [3] o ~The algoritmums so developed may be described 

as procedures that make maximal use of minimal cost routes, Each assigns 
quotas only to the least expensive feasible routes until the assignment 
is complete. The resulting assignment is one which gives minimum cost. 
It is not necessarily unique. These algorithas are also discussed more 
fully im Appendix A and a problem is solved using each 1m turn, 

3, The Transhipment Technique. 

An interesting variation on the m x n transportation problem (m 
sources, m destinations) is the transhipment problem described by Alex 
Orden [9] In the original transportation problem each point acted as 
a shipper only or as a receiver omly, and the routes were considered to 
be direct from each shipper to each receiver. Im his extension of the 
problem, Orden has permitted shipments to go wia amy sequence of points 
rather than being restricted to direct connections from one of the origins 
to one of the destinations. His method of solution involves converting 
the extended problem to the original problem by a rather simple device. 
Each point is treated as a pair of points, ome acting as a shipper and 
one as a receiver, The unit cost of shipment from a point considered 
as a shipper to the samé point considered as a receiver is set equal to 
zero, All other of (i, j =1l,..0., m+n, ity) are assumed known 
and greater than zero. For computational purposes large stockpiles are 
assumed to exist at each shipping point, The excess of stockpiles over 
actual amounts shipped appear as shipments from a point to itself at 
zero cost and aré removed from the final result. With this mechanism, 


the mxn problem with transhipment is converted to an M x M problem 





without ttanstipmenmt «qa Hirchcack problem), where M=m+n. “The 
» é 
problem 15 Solwed tw this fore ard, im the non-degenerate case, 2M 


- 1 routes have positive quotas assigned. It nappens chat in minimum - 


cost solutions ail «x, are positive, but they are removed leaving M - 


t 


Lk 


lemtn-l actyval shipments, 
4, Integer Solutions. 
It is worth noting here that if the Hitchccck problem is posed in 


integers (i.¢. each a, and b iS a positive integer) then there is 


i j 
at least one minimizing solution im which each « is either a positive 


i] 


imteger or zero L4 ] o in this case the whole problem can be solved in 
integers by introducing triai, feasible solutions consisting of integers 
and then changing these only by integral amounts. Very recently, Gomory 

[13] has developed a “method of integer forms’ tor solving a general 
linear programming problem with the additional constraint that the 
variables in the solution be integers, Dantzig p14) has included this 
procedure in the term “discrete programming” and forsees in its further 
development promise of solving “all kinds of combinatorial, non-linear, 
non-convex problems..,.areas where classical mathematics has been weakest". 
5. The Capacitated Problem. 

In the classic Hitchcock problem the upper bound on the amount of 


commodity assigned ta amy route is determined by the amounts available at 


ft 


i 


The routes themselves are not regarded as Limited in their capacities. 


the sources, a, $, and the amounts required at the destinations, Ps 


One can imagine practical situations wherein the actual capacities of the 


routes may be limiting factors. For example, suppose it is desired to 


*See Appendix A for a discussion of degeneracy, 





transporc cértain known quantities of commodity trom several sources to 


7 y -*, ae ae, p? uae * - = re 46 « - a ¢ 
everal destinations each dey, or wonth and we rave available certain 


(ha 


railway rovtes with specified capacities; then, the amount of the com- 
modity shipped om route i-43, in the specified tims interval, cannot ex- 
ceed the capacity of that route. 

Our introduction of capacities has trougnt time into the problem, 
at least to a limited degree, The natural] specification of capacity is 
arate, Consider @ problem similar to cme Hitchcock problem, except sup- 
pose the a's répresent quantities to ba shipped im a certain time in- 
terval, say a month, and che Pe are amounts required at the destina- 
tions per month. Further, let us specify che capacity of each route i-j 
im umits cf commodity per month, We have then the ¢o - called capacitated 
Hitchcock problem. This preblem iS Simisar to the uncapcitated problem, 
except thac we envision it as being repeatedly applicable within time 
sintervais of specified length. Basically, it is static in that we have 
introduced no truly dynamic variables such as instantaneous flow rates 
or transic times. And we have suppoged the given COMSLANES- 4, S, b's, 
and the capacities- do not change with time, 


We state the problem formally. 


a) 


Given. 
(a) Two sets of constants 
a i = Lyocs ym 
j j 22 eee 
Such that = a, 


Ne 

er 
Ci. 

bi 

om 


(b) A matrix C « Ee aan! 


‘¢) another matrix oD - | 


Le 





where “4 represents the capacity of routs ij. and 
ag represents the cost of transporting & unit 
amount from source i to destination j. 
Find X = [ox J 
Lj 
such that 
(1) 
(2) 
p 
(3) Of x., = c,, 
(4) S oa x, ,d. = minimum 

Whereas the uncapacitated problem always has a feasible solution C4) 2 
the capacitated problem may not. The problem is feasible if, and only if, 
the maximal flow of the associated network is at least equal to the total 
amount to be shipped, K. A discussion of networks and their maximal 
flows will follow in Chapters II and IIT. 

Ford and Fulkerson have also provided a method for solving the capa- 
citated problem in their paper, "A Primal - Dual Algorithm for the Capa- 
citated Hitchcock Problem” (is . As their citle implies, these authors 
define a dual problem involving dual variables ee fo } and 0 44 
These are used to distinguish certain routes, and, in terms of them, a 
restricted primal problem is solved. From this solution, new dual vari- 
ables are defined and the restricted primal problem is again solved, this 
time using the new dual variables. Im each solution of the restricted 
primal problem, a set of positive Ae ie determined. When the sum of 
these ee is equal to the total amcunt to be shipped, kK, the 7 SO 
determined provide a solution to the original problem, 


During each iteration, in determining the solution to the restricted 


oO 





primal problem, the authors call upon a special tlow algorithm developed 
by them in 4a previous paper [26] . The mechanics of tie solution are 
Similar to these described for the wncapacitated problem and will not 

be described in further detail here. 

In a subsequent paper (17] » Fulkerson showed that an upper bound 
on the number of iterations or "labelings” required to solve the trans- 
portation problem could be determined. These upper bounds were: (a) 
for an mn xn optimal assignment problem, n(m+l); (b) for the tran- 
shipment problem, K(n+1); (c) for a Hitchcock - Koopmans problem 
K [. min(m,n)+1_] : 

The literature contains numerous other examples of applications of 
linear and non-linear programming techniques to transportation problems. 
Several such papers are cited in the bibliography without special reference 
to them in the text. These are intended ts be representative, rather than 


exhaustive; to attempt an exhaustive Listing would be presumptuous, indeed, 


LO 





10, 


li. 
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CHAPTER If 
TRANSPORTATION NETWORK'S AND GRAPH ‘THEURY 
1. Dauscussion. 

Tn the study of transportation systems it ts often dasirsble ta 
utilize network diagrams, The usual approach is to let the nodes of the 
network represent origins, destinations, or junction points. and to con- 
mect these nodes with lines representing the available routes. The des- 
cription, representation, and use of such diagrams has been greacly facili- 
tated by the use of the mathematical theory of lintar graphs An early 
example of the use of graph theory in connection with transportation 
problems was in a paper by T. C. Koopmans and S. Reiter [1] published 
in 1951, Much of the work had besn done earlier during World War ITI 
by the former duthor. 

The time - lag in the application of graph theory to transportation 
problems is emphasized by quoting from the above mentioned paper. 

"The cultrual lag of economic thought in the application of mathe- 

matical methods is strikingly illustrated by the fact that Linear 

graphs are making their entrance into transportarion theory just 
about a century after they were first studied iam relation to elec- 
trical networks, although organized transportation systems are 

much older than the study wf electricity.” 

Perhaps some of this delay is due to the fact that, uncil wery recent- 
ly, graph theory has been somewhat neglected by mathematicians themselves. 
Few articles on the subject appear in mathematical journals. 

Until 1958, there was only one book devoted entirely to graph theory 
- D. Konig's Theorie der Endlichen wid Unendlichen Graphen [2 } » first 
published in Leipzig in 1936 and reprinted by Chelsea in New York im 1950, 


The second book on the subject, by Claude Berze (3) , was vublished im 


13 











France in 1952, 

Professor F. Harary, of the University of Michigan, has written 
quite extensively on graph theory and, with R. Z. Norman, is writing a 
book due to appear in 1960. [5 | 

Because of its wide applicability and comparative simolicity, the 
two-terminal . network is used in many transportation problems, Normal- 
ly, flow in such a network is unidirectional from a source to a sink. 

This theory has been especially useful in determining the maximal 
flow which a network will sustain. Examples of its use in this con- 
nection are included in Chapter [II. 

The similarities between communications and transportation net- 
work problems have been recognized and pointed out by R. E. Kalaba and 
M. L. Juncosa C4] » among others. This similarity is jot 14wigee to 
their common tie with graph theory but extends also tto queuing theory, 
Simulation devices, linear programming, dynamic programming. and Boolean 
algebra. 

2. Representations. 

We shall exclude from our discussion graphs containing arcs which 
join individual nodes to themselves. 

A natural representation of a network is a diagram in which points 
represent nodes, and lines represent the routes between them. Elsewhere, 
(Appendix C) we have defined a network to be a (connected) graph G, to- 
gether with the capacities of its individual arcs, We shall illustrate 


here 4 oneé-origin, ones terminal planar network with source a and sink b. 


1, For definitions of this and other technicai terms, see Appendix C, 





The graph 1s net oriented. 





5 ° ee Hh, 
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This network covld ajJso be represented by &@ symmetric matrix 


I =(2..| in which  & . is the ctpaciry of arc uj if the are 13 
13 1j 


is included in the graph, and is zéro otherwise. For our example the 
applicable matrix would be che tne velow, Except on the main diagonal, 


heros are not written but are implied by a blank Space, 
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In Chapters ILI and IV we shall wse the lanzuage or tie theory of 


graphs, For representations of the grap’. sich we go2cy dlaprams such 


\ 
z= 
r 
Cc 
4 


as figure 2-1 rather than the matrix representation, firure #72, will 
be used, 

It is desirable to emphasize that we have pointed cut only one field 
for the application of graph theory. In addition to their usefulliness in 
the study of transportation, communications, and electrical networks, 
graphs have been used in many cther fields. They are useful in describ- 
ing Markov chains, for example. As another example Professor Harary 
has made use of them in the study of the behavicral sciences C6] q C7] : 
Much of his work in graph theory has been done in conjunction with the 
Research Center for Group Dynamics, University of Michigan, and the 


Social Seilence Research, Bell Telenhone Laboratories, 
> a 


16 
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CHAPTER Iit 
MAXIMAL NETWORK FLOW 
1, General 

As we mentioned in Chapter II, transportation systems are often re- 
presented by metwork models. The nodes may represent simple transfer 
points such as the pick up point for the hook in a cargo-handling opera- 
tion; or, they may represent highly complex entities such as a seaport 
or a complete railway operating division. What the nodes represent de- 
pends upon the model under consideration and the degree of generality 
or detail under which it is being studied. These nodes are joined by 
arcs which may have an associated traffic-carrying capacity, This 
collection of points and arcs with their capacities will be referred 
to aS a network - (see Appendix C). 

One question of interest is that of determining the maximal 
steady-state flow through such a capacitated network. We shall des- 
cribe in this chapter two procedures for determining maximal flow, one 
a flooding technique due to Alexander W, Boldyreff [1] and the other 
the minimal cut theorem dué to L. R. Ford, Jr., amd BD. R. Fulkerson (a 
2. The Flooding Technique 

Boldyreff applied his technique directly to traffic through a rail- 
road network, and was careful to claim a6 generality for it. He character- 
ized it as empirical in nature, amd left as an open question the possibil- 
ity of wider use of the flooding technique. It appears feasible to apply 
it in any situation where a network model will fairly describe the actual 
situation. 


The method employed can be applied to networks with any number of 


! 
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origins and terminals, however, for ciaricy the one-origin, one-ter- 
minal case is described. 

Certain subsets of points are distimguished a5 the origins, the 
terminals, and junction points (see Appendix C). We will discuss flow 
in a4 railway network assuming capacities in trains per day. By ‘steady- 
state"’ wa mean that the flow rate through each arc, in trains per day, 
is constant from day to day. The steady-state is further clarified in 
these initial assumptions: 

1. At the origins omly loaded trains leave and only empty trains 

arrive, 

2. At the terminals only loaded trains arrive, and only empty 
trains leave. 

3. At each junction point the number of ivaded (empty) trains 
arriving is equal to the number of ioaded (empty) trains 
leaving each day. 

4, The number of loaded trains leaving all the origins is equal 
to the number arriving at all the terminals each day. 

a The mumbexr of empty trains leaving ali the terminals is 


equal co the mumber arriving at ali the origins each day. 


The last two statements are deducible trom the first three and are stat- 
ed explicitly only for emphasis, 

The problem is to find a methed of assigning the steady-state 
flow of traffic to each are of the network, mot exceeding the capacity 
of the arc, which will maximize the total flow of loaded (eampty) trains 


from the origins (terminals) to the terminals (origins), and satisfy 
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the five assumptions =ttated above 

In the problem as stated, the *teady-state consises of simultaneous 
movemenr through the network of loaded trains i: rum the origins to the 
terminals and empty traims in the opposite direction, These two move- 
ments must necessarily be at the game rate. By imagining a network 
with are capacities just one-half of those cf che original network, we 
can look only at the flow of loaded trains from origins to terminals and 
reduce the problem to the unidirectional case. By unidirectional flow we 
mean that the flow in any one arc is im one direction only. In this dis- 
cussion the flow considered will be generally directed from the origins 
to the terminals. 

For this simplified problem the steady-state conditions are con- 
tained in the following statements: 

1, At the origins all trains leave, nome arrive. 

2. At the terminals, ali crains arrive, mone deparcc. 

3, At each Junction point the number of trdins coming in is 


equal to the mumber going out. 


+ 
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4, The number of trains leavin l origins must equal the number 


r pa’ 
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arriving at all terminals. 
5. The traffic flow through each arc cannot exceed the capacity 
of that arc, 
By defining two networks to bea equivalent if cne maximal flows chrough 
them are equal, we can readily reduce # network with many origins and 
many terminals to an equivalent network with a single origin and a 
Single terminal. 


We shali describe the flooding technique utilizing a sirele-origin, 





single-terminal metwork. The procedure is as follows: 

Starting at the origin, assign sufficient traffic flows to all arcs 
leaving the origin to saturate them. This gives the maximum number of 
trains arriving at junction points one arc removed from the the origin. 

View this set of junction points as new origins and, starting with 
the one subject to the greatest capacity constraint, schedule trains 
whenever possible in the following order: 

Ll. 'Forward’ - to new junction points through the outgoing arcs 

2% "Laterally’ - to other points of the set. 

34 "Bottlenecked’ - if trains are left over after steps 1 and 2; 

i.e. if all outgoing and lateral arcs are saturated. 

At such a set of junction points there may be arbitrary decisions 
regarding which trains to ‘bottleneck’. The guiding principle is to 
move forward as many trains as possible, and to maintain the greatest 
flexibility for the remaining network. 

Continue the above procedure until the scheduled flow covers the 
complete network and reaches the terminal. 

Eliminate bottlenecks by returning all excess trains to the origin. 
The validity of the solution can be checked by inspection. If a maximal 
flow has been achieved, there will be no continuous umsaturated path 
extending from the sna ibe the terminal. (That is there will be no 
chain Perec t inguokivin and terminal that does not contain at least one 
saturated arc), If this criterion is mot met, the inspection will have 
revealed the unsaturated chains. Flow in these chains can be increased, 
giving maximal flow. As will be seen, the observation concerning satur- 


ated arcs ig a point of similarity between this procedure and the 


2i 





misimg. cul procedure 25 *@ discussed, 

In this paper, woldyretf also gives some procedures for simplifying 
complicated networks. However, as he points out, the networks to which 
these procedures can be applied are rather more the exception than the 
rule, and the applications are tedious. The author recommends a straight- 
forward applicatiom of the tlooding technique, 

A sample problem is solved ty this procedure im Appendix B. The same 
problem is also solved using Ford and Pulkexgon’s minimal cut theorem. 

3. The Minimal Cut Theorem. 

This method was applied to networks, um general, not necessarily 
rail networks, Although we will not give here &@ proof of the minimal 
cut SheokEm’, we will meed some defimitions in order to discuss its ap- 
plications, See Appendix C. 

A graph, G, is @ fimite, ome-dimen¢gional complex, composed of 

vertices a,b,¢,...,¢ and arcs A (als), B Cla ee d (ce). 

An arc lan (ab) joins its end vertices a,b it paeses through no other 
vertices of G, and intersects other arce only im vertices. 
A chain is 4 set of distimet arcs of G which can be arranged as OX (ab), 
3 ‘oe od igh) where the wertices @.6.¢,...,h are distince. A 
chain does not Intersect itself, it joins its end vertices a and h. Each 
arc in G has jecéc Tamme wer it a positive number called its capacity. 
The graph G, together with the capacities or its individual arcs, is 
called 4 network. Two vertices of G aré distinguished a, the source and 


b the sink. A chain flow from a to b is & couple (C,k) composed of a 


l. For @ proof cf the thorem, see tne second article lizted in the 
bibliography at the end of this chapter. 
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chain joining a and b and a non-negative wumber k representing the flow 
along C form sowrte to sink. A fiow in a network is a collection of 
chain flows which has the property that tne sum of the numbers of all 
chain flows that contain any arc is no greater than the capacity of that 
arc, If equality holds, the arc is said to be saturated by the flow, 
A chain is saturated with respect to a flow if it contains a saturated 
arc, The value of a flow is the sum of all the chain flows which compose 
it. A disconnecting set D, is a set of ares which has the property that 
every chain joining a and b meets the collection, A disconnecting set, 
no proper subset of which is disconnecting, is called 4 cut. 
The value of a disconnecting sat v(D) is the sum of the capacities 
of its individual members, Thus a disconnecting set of minimal 
value i¢ automatically a cut. 
We state the minimum cut cheorem. 


Theorem 1. The maximal flow valué obtainable in 4 network 1s the 
minimum cf wD) tcaken over all disconnecting sets |. 


This theorem is so intuitively appediing that it hardly ssems to 
require proof at ail. Wowever, 4&8 examination of reference [2] will 
Show, the proof is quite involved, ani for certain networks (for example 
see C3} ) it does not apply. The inapplicability to networks with 
Several sources and corresponding Sinks wae illustrated also by Ford and 
Fulaerson (2 ] : 


The actual computation procedure for determining the values of the 


on 
FF 


minimal cut is based om a corollary tc the minimal cut theorem and an 


additional theorem which will be stated without proc. 


Corollary: Let A bé & collection of Sres of a network N 
which meets «ach cut of N in gust cme arc. Tk WW’ is a net- 
work obtained trom N by adding « “re capatity of @ach 
azc of A, then cap (N') = cap oMe + & 
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Theorem 2, Lf the graph of CC, "“ogether witn arc ab, is a 
planar grapn (4) . chere exists a chain joining a and b which 
meets each cut of N precisely since. 

Let T be the chain joining a and b which is tcpmost in NH. T has the 
property of theorem 2, Impose as large a chain flow as possible (T,k) on 
this chain, thereby saturating one of its srces,. 

By the corollary, subtract k from the capacity of each arc of T. 
Remove the previcusly saturated arc whose capaciity is now zero. Continue 
this procedure. Eventually the graph disconmects and a maximal flow has 
been established. 

For an illustrative example, see Appendix 5. 

In problems of sea-transport, the most restrictive capacity limitations 
are quite often the cargo-handling capacities of ports ar which ships must 
load and/or unload, rather tham on sea rowtes themselves, Tnis is often 
true «ven when the number of ships availatle is limited. &t first, it 
would appear that this is a preblem wherein the capacity restrictions are 
on the nodes of a graph rather than its arcs. We Sugeest, howaver, that 
this can be reduced to the previous problem by employing a device similar 
to that used by Orden. See Chapter I. Treat each node representing a 
port of limited capacity as both a@ receiver and a shipper and replace the 
single node by two nodes joined by am arc showing the capacity. We have 
merely introduced a slight modification to the model to shaw more detail. 

The applicable units of capacity might be troublesome in a ship- 
ping network due to the wide variety of ships with which one mignt deal. 
There is 4 great temptation to dismiss this difficulty by introducing 
some sort of 4a “standard” or “notional” ship, and for a first approxi- 


mation this may suffice, However, some ports may well be able te handle, 
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for example, 3 notional ships per day provided this capacity actually 
came in three ships near our “standard” size, and be fineapable of even 
docking one very large ship of equivalent capacity. 

Thus our model begins to break down anc we need to introduce a more 
realistic one capable of handling the real situation. We would emphasize 
chat the models with which we have dealt to this point are mostly systems 
models dealing with large-scale inputs and gross effects. They would re- 
quire inputs representing aggregation of many lesser inputs, and these 
must often be combined by skillful application of lessons provided only 
by experience. Simple arithmetic summations rarely suffice. The inter- 
pretation and application of the results prowided by the model is equal- 


ly difficult, and requires the same careful judgment. 
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MINIMUM TRANSPORT UNITS TO MATNTAIN A FIRED SCHEDULE 
1. General Discussion. 

We shall consider in this chapter a problem which 18, im 4 Sense, 
complimentary to the problem of maximal flow previously discussed, We 
shall be concerned with the determination of the minimum riumber of trans- 
port units necessary to maintain a fixed schedule. 

In order to make precise the notion of a schedule, we shall again 
use the language of graph theory. We envision a conmected graph in which 
a subset of nodes is distinguished as terminals; others represent junc- 
tion points, Sy a run we mean the traversing of a certain path from one 
terminal to another by a single transport unit. Runs originate idepart) 
and terminate (arrive) at certain specified times within an interval, 
such as a day or a week, These departures and arrivals are repeated in- 
definicely in past and succeeding intervals of the same duratiscn, We 
choose a typical time interwal of length 77. commencing at some arbi- 
trary time t and ending at time t) = t) +77. & schedule is the 
specification of the time of each departure and arrival during the time 
interval 77 , tor each terminal in the system. It will be noted that the 
duration of a rum may be less than, equal to, ox preater than ouzr typical 
interval, 

Perhaps the simplest feasible solution cf a schedule problem is to 


assign each unit to repeated round-trips b 


tty 


tween two terminals, This may 
involve the movement of @ great many empty or "dead-head” umits, and is 
rarely optimum, if other utilizatiom of units is feasible. As can be 
readily seen, the problem of minimizing transport wnits to weet a pre- 


scribed 


& 


nd 


chedule reduces to the ¢roblem vf routing empties ri) 
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We have found several procedures for determing minimum units for 
fixed-schedule operations, L2] ; (3 ] | 1 (s} . in this chapter we 
shall describe, as a representative example, an algcrithm by T. E. Bartlett 

[3] built around a railway schedule. 
2. Algorithm for Minimum Transport Units. 
Assumptions 

(1) Schedules under consideration refer to rums of transport units 
between terminals. The runs originate or terminate, at specified times, 
within an interval, of specified duration, and are repeated indefinitely 
im past and succeeding intervals of the same length. 

(2) Each terminal has at least one route connecting it with one or 
more other terminals. Any terminal may be reached from any other in the 
system, The system may be represented by a connected graph. 

(3) Each terminal has as many “departing"’ runs as “incoming™ runs. 
This may require transfer of empty or "dead-head" units. 

(4) For gach run there is a standard or minimal type of transport 
unit which is employed. Units may be employed on any run for which speci- 
fications are met, 

The above assumptions are quite reasonable and are usually satis- 
fied by transportation systems operating on schedules, 

Additional Assumptions 

For the particular algorithm presented, certain adaLttional, restric- 
tive assumptions were made. These were 

(1) Requirements at any terminal for transport units to make up out- 
going runs must be obtained from regularly scheduled incoming runs at that 


terminal, 
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(2) All transport units are of a single type. 

(3) All schedules are met precisely. 

(4) No delay beyond the time differential of scheduled arrivals 
and departures is required for interchange of transport units from one 
run to another. 
Discussion 

By the hypothesis concerning the schedule a umit is either on a run 
or it is “idle” at some terminal awaiting transfer to 4 run. A schedule 
which is maintained generates a total running time plus idie time which is 
equal to the number of units employed multiplied by the leneth of a typical 
period. Since the total running time is fixed by the schedule, the number 
of units employed will be a minimum if, and only if, the total idle time is 
at a minimum. Therefore, what is required 14 4 pairing of arrivals and 
departures at each terminal which will minimize toral idle time. 

Three séparate cases are considered, 


Case I, All departures in 4 typical period at a terminal are 
later than all arrivals during this period, 


Case If. Some departure is earlier than some arrivals, but all 
arrivals and departures within the period can be paired 
in at least one way so that each departure is later than 
its paired arrival, 


Case IIL. Up to some time during the period, more departures chan 
arrivals have been scheduled. 


In the first case any pairing of departures and arrivals generates 
the same total idle time. In the second case, pairing of the "first in" 
with the "first out” is both feasible and yields minimum idle time. 

Case LIT presents a special problem. Since there are, up to 4 point, 


more departures than arrivals, there is no feasible pairingof all departures 


with arrivals within ons period, It 15 necessary to pair one er more 
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departures with arrivais from a preceeding period and some arrival (3s) 
within the period are left over for paizing with departure (sjin a sub- 
sequent period. By extending the sequence both forward ard backward 
sufficiently to pair all departures with arrivals, the pairing problem 
reduces to that of Case ITI, and the first-in-firec-c.ut procedure gives 
minimum idle time. However, it is mecessary co adept a consistent assoc- 
tation of idle time to a particular period in order not to miscownt the 
total time. The procedure adopted was te associate with the current per- 
iod that idle time generated by pairing runs of the previous period to 
runs in the current period as well as the tile time generated by pair- 
ing within the period, 

Idle time developed in the pairings from the current period to the 
next period is not charged to the current period, but to the next one. 
Results 

From the above considerations 4 toral, remning plus idle, time is 
calculated, and ig shown to be. 

te WiK+*Z BD 
where 


Iv = length of period 


of = total rums en roure at the be erinning dor 
end) of a period tor 2ll terminals 

ef = maximum value cf cumulative departures Less 
cumulative arrivals at terminal j during a 


pericd, 
The minimum number of transport wnits required is then obtained by 
dividing T by the length of the period 77” , and is giwen as: 


P 


U = minimum number of transport units = 


Ss NA 
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In concluding his paper, Bartlett pointed cut some areas of schedul- 
ing problems under study, Thes¢ are listed kelow, Some wf them, if can 
be seen, are aimed at removing the restrictive sssumpticons made earlier, 
Items listed were; 


a. Determination of routings to parmit economic maintenance 
schedules, 


b. Determination of routings to optimize the numbez of units 
used for the situation in which several classes of equipment 
may be employed on various routes, 

c. Determination cf routings to minimize the number of units 
used and maintain schedules having probabilistic arrival 
times, | 

d, Determination of optimum schedules based om cost of trans- 
port umits, operating costs, and demand distributions, 


Eo Analysis of the possible application of the concepts in 
the field of machine scheduling. 


The actual problem of allocating units to 4 schedule im a manner to 
preserve minimality and, at the same time, meet certain maintenance require- 
ments, was discussed in a subsequent paper by T. E. Eartlett and a. Charnes 

[5] » Their method of solution involves translating the schedule into an 
oriented or directed graph, and working with tne associated incidence 
matrix. A node in the graph considered represents a specific time at a 
specific terminal, The time chosen is che time of a local maximum in the 
"idle equipment inventory" at the terminal under consideration, There may 
be more than one "node" for each terminals there will be at least ome. An 
arriving rum prior to a node time, which may be paired with a departing run 
after the node time, without violating minimality of equipment, is regarded 
as a link (arc) positively incident on that node, A departing run after 
a node time which can be paired with an arriving run prior to the node 


2 


time is viewed as a link negatively imcident on the ncde, This representa- 


tion gives seme icea ct the versatility of graph theory in the study of 
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transportation probiems. 

With the incidence matrix constructed a5 described above, each 
coiumn corresponds to one of the runs of the system, and each row corres- 
ponds to one of the nedes. Ome can then trace out a cycle, or a "pool", 
in the following manner. Connect 4 minus one with a plus ome in the 
same row, and repeat, connecting pius ones and minus ones alternately 
within the same column and then the same row, When the first or start- 
ing emtry is connected the second time, a cycle is formed. The pairing 
of runs in this manner is both feasible and maintains minimality of 


equipment. 
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CHAPTER V 
MISCELLANEOUS TECHNIQUES 

Up to this point we have dealt with a rather restricted group of 
transportation problems and models. For the most part the models hinve 
been of the network type and descriptive of an entire transportation 
system In this chapter we shall give examples of the application of 
other techniques te rroblems in transportation. Specifically we shall 
cite examples using queuing theory, the Monte-Carlo method, and engineer- 
ing analysis. 

Among the earliest models employed in the study of transportation 
systems was the queue, or waiting-line, model. Examples of queues in 
transportation systems might include automobiles at a stop sign or a toll 
station, freight cars in a freight yard, airliners circling an airport 
awaiting clearance to land, or ships in harbor awaiting unloading. In 
a different but related field this model has been applied to telephone 
calls at switchboards and messages in communications relay stations. 

Waiting line theory is strongly related to probability theory, 

Since both the arrival times of ''customers" and the service time for an 
individual customer are stochastic variables. This is true to a large 
extent even though arrivals might be "scheduled". The usSual practice is 
to study the so-called steady state; that is, one assumes the queue has 
been in operation long enough for transient effects to have diminished to 
anegligiple point. One usually knows, or can Soon estimate with accep- 
table accuracy, the mean arriyal rate of customers and the mean service 
time The specification of the distribution of these random variables 


is more difficult. However, the most widely used assumption is that the 
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number of arrivals within a specified interval bas 4 Poisson distribytion 
and that service times have a negative exponential distributicn. Witn 
chese assumptions and inputs, the mcedel can be analyzed to predict medn 
holding time - waiting time plus service time -, mean queve length, and 
the probability 4 customer will have to wait. 

& good example of the application of queuing thecery to a transporta- 
tion problem is contained in an article by Leslie C. Edie cil of The Porc 
of New York Authority, in which he discussed traffic delays at toll booths 
of several tunnels serving the city of New York. 

In many problems involving queues, sclutions in closed form are 
Gifficule, t£ not impossible, to obtain. Such problems are sometimes amen- 
able co 4 certain amount of analysis by the Monte-Carlo method. In using 
this m@thod, one develops a model of the system he wishes to study and 
then aimuiates arrivals from time to time. The arrival times are assign- 
ed at random in accordance with previcusly observed arrival statistics. 
Nirh che model the investigator can determine the effect of varying certain 
—atamercors in the system. He can, for instance, introduce an additional] 
Service channel or allow priorities. With the simulation technique, the 
researcher can ger data corresponding to months of actual operation quick- 
ly and cheaply. 

Roger R. Crane, Frank B. Brown, and R. O. Blanchard (2] have reported 
on the application of this method to a railroad classification yard. These 
authors developed a model by the method described above and, using simpli- 
fied statistical assumptions, isolated areas for more careful study. [In 
the two areas reported they determined the effect of substituting parallel 


yYocessing for series processing and of changing the "rules'' regarding the 
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lem gf awitch poe .res considerable savings were reaiized in beth im- 
Scan ‘PoyceSsing in this context consisted of inspectios ard cimsesi-+ 
fication. 

For our last example we nave chosen the Cargo-Handling kesearch 
Froject in the Department of Engineering, University of Celifornia ac 
los Angeies. This research 18 sponsored by the Office of Naval Research 
and the Maritime Administration, Though by no means limited thnereta, 
this preject has made use of engineering analysis, statistical methods, 
and the Monte-Carlo technique. The project has been reporred in derail in 


he &@ 


=¥ 


a series cf numbered reports, (references (a ] ; [5] s (5) - C7] Ne 


Ch 


are not generally available. R. R. O'Neill reported on the project an 
discussed the analysis and Monte-Carlo simulation of cargo bhandline bde- 
fore the Purdue University First Transportation Research Symposium, in 
February 1957, This presentation, along with others given at che synpo- 
sium, was published in Naval Research Logistics Quarterly [3] im September 
eos’. 

fhe transfer of cargo between ship and shore has long been recognized 
as the major bottleneck in the operations of the shipping industry. The 
basic objeccive of the research at the University of California at Los 
Angeles has been tne effecting of improvements to the cargo-handling 
SyYSter. 

For the analysis of cargo handling, the project faced the problem 
ef choosing a level of description which would yieid to analysis and yet 
lead to useful results. The level chosen considered the cargo-handling 
system as 4 series of transporting links in which each link was the sum 


of aii the movements required of a carrier in transferring a load cf 
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cargo from one point to another. In a land-to-water sequence the links 
studied were: (a) the shed link, (b) the wharf link, ic) the hozk 
link, (d) the hold intermediate link, and (e) the hold link 

With the cargo-handling system idealized in the foregoing manner. 
certain hypetheses were made concerning the interrelations existing 
among the various links. Two examples cf such hypotheses are. 

Hypothesis I: [If there is a group of transporting links connected 

in series, then there is one link which has the least mean induced 

delay and is the slowest or controlling link. 

iypathesis IL: The relationship between the relevant factors which 

affect the mean activity time of a link can be formulated se as to 

make prediction possible. 

Both of these hypotheses were used as the basis for a fieid study con- 
ducted in the summer of 1953, primarily in the ports of Los Angeles and San 
Francisco, Even though a random, work-sampling technique was used te obtain 
data, there were 41,218 individual observations recorded, This gives some 
idea of the magnitude of data-gathering problems. 

it is noted that at the gross level a port terminal can be thougnt 
of as a node with a land carrier as the link on one side and a ship as 
the link on the other side. The approach used in cargo-handling study 
iliustrates the meaning of a remark made earlier concerning the necessity 


of building models in a degree of detail consistent with the problem to be 


Studied. 
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APPENDIX A 
ALGORITHMS FOR THE HITCHCOCK DISTRIBUTION PROBLEM 
1, The Dantzig Algorithm. 

In this appendix, we shall describe the essential elements of Dant- 
zig’'s algorithm for solving an m x n transportation problem based on his 
Simplex method. See Chapter I. We shall then work out a simple problem 
uSing this method, We shall make similar presentations of the algorithms 
of Ford and Fulkerson, and Munkres, working out the same problem in each 
instance, The problem is 4 x 5 (four sources, five destinations) and the 
numbers are quite simple. We hoped to choose one which could be followed 
easily and yet was sufficiently difficult to illustrate the methods. The 
problem is stated below: 


There are four sources with amounts to be shipped as follows: 


a, = 10 
a, = 60 
a, = 40 
ea 90 


There are five destinations with amounts required as follows: 


by =i) 
do = 40 
by = 70 
Dy = 30 
be = 40 


We note Za a, = os b, = 00 


a7 





The cost of transporting a unit amount of commodity from source i to 


destination j is given as the ij entry in the following matriy: 





Our prceblem is to find: 


X = Cx, J such that 


* 4,20 eee 3 404 2125384 5 iol 
S my Ay eri 253.4 ey 
3 

i J = 1,2,3,4,5 1.3 


aoe x. = Minimum 
Before commencing the algorithm, we shall state, without proof, some 
useful theorems relating to the problem. For proofs the reader is referred 
to Dantzig fi] and Gass [2]. 
Theorem 1. The transportation problem has a feasible solution. 


Theorem 2, A solution of at most m+n - L-positive oe 


exists. 


Theorem 3, If the a. and b. are all non-negative integers, 
then every basic fdasible solution has integral values. 


Theorem 4, A finite feasible solution always exists. 
For 4 more general case with m sources and n destinations, if we write 


out the equations 1.2 and 1.3 we have m + n equations in mn unknowns. Eut, 
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Sivce we have specitied > a. 


and the system may be reduced tom+n- 1 singarty indpendent squatiruns. 


A b@sic feasible solution then has m+n - 1 positive ie fh tise 


feasible solution is found with less than m +m - 1 positive X, ‘®, deg- 


anuracy i¢ indicated. The difficulty can be averted by perturbing the 


phi ter an the following manner, add a small positive amownt € to 
swe) 2 leave ali b_"s unaltered except for b , and add n€ to b 


od 
u 


inle preced: re need not be employed until degeneracy is indicated, Solve 


7 


the frrbien for a minimum, basic feasible sulution, using the € 5 if 


Atter the sclution is obtained allow each € to go to zer>. The sclu- 
Cien 18 still valid. In our problem, since «ach ai and OF 1552 posi- 
tive imteger, we are assured of a solution involving cnly positive in- 


c 


ip 


VES, 

To €xplain the computational procedure we shall work out cur problem 
Step by step. As the initial step, we write down the matrix cf direct 
cost cosifticients and a blank asSignment matrix. We shall habitually 


rite the column of a,'8 to the right of the assignment matrix, end a row 


2 


cf b,'s below the matrix. This procedure will serve to keep the pro- 
3 ° 


blem before us, 
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bee bye fee gonence Winimurn [ 40,60 5 10 | ~ GO Sontanuing in this 


Pay. we Ghtsain the following basic, feasibie feolution. 


An 
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In tnis example we have mtn - 1 = 8 positive X'S lf afcer 


Cone 


having madé one cf our assignments xy the amount remaining in the cnt- 
umn j and that remaining in the row i had vanished simuitaneously, we 

would mave had a degenerate solution, except for the case i 2m and j = 
m. In case of degeneracy, inclusion of the € perturbation, as described 


above. will remove it; for then no partial sum of the a's anid ae can 


equal any a OF [ogee 
J 


For comparison purposes we note here the cost of our first basic 
feasible solution. 


Cy * 1014) + 10¢2) + 4093) + 10(5) + 4043) + 20¢4) + 301) + 40(2) = 540 


Next we compute, for the basic feasible solution just obtained, its 
indirect cost matrix. Each element ak of che indirect cost matrix wiil 
be the sum of an indirect cost us associated with the source (row) and 


an indirect cost ve associated with the destination (column). Thus 


o. = wu, + ¥, for all i and j. 
ij L 3 


For those pairs i j employed in the trial solution 


Comm Ce. 

ij ij 
From this latter equation we can get some of the entries for the indirect 
- cost matrix. In our example we get the incomplete matrix below. 


Indirect - Cost Matrix (incomplete) 

















we cdnan@e ifr ou & 4 andy v, = GO, arblQwer.iy) 287 Goemerg ete ve- 


} a . . Cm - - -” - . = 
MAL UY eG > @ei 7 3s f£=am tae we ataby 
- ii 


c Pee 4 for ali 1 ana § 
~3 “ j 
Doing 32 w+ obcaige ch indirect cott mAtrin Below. Tce ; . 


hence by the equation U, + 4) = Cy, > as Wy mast equal 2. Then ¥, ™ i 


and ¢ om, until ail uw, and ae are determined The remaining Bn are 
a 
comover i by using the relation, c¢c,, = wu + ¥Y,. 
ij L i 


imdirect Cost Matrix 





In order ta facilitate comparison of corresponding entries in the 


diz 


(Ts 


et and indirect cost matrices. we shall write them together as one 
matrix, called che cost matrix. In each cell of the cost matrix we 


sm@il write ¢ in the upper &-1f and ¢.. in the lower hali, In our 


ty 3 
sample problem this reiJits .> L116 foilowing cost matrix for step fb. 


1d 














neriec above thre diarwal 


Batries below che diagenal 


The acditircnai column ™" CFF 


The accirctional row below ts 


cost Matrix 
(step 1) 


are Gee i 
J 


i 
p< 
ie 


a ee Y 


—~ame 





——— 





We now leok at those cells in which ae 7” Cc Tf there are none 


ij’ 


the basic feasible selution we have chosen is a minimim cost sTantiga and 


our problem is solved. If C14? C5; for some 1 and 4, we wael te 
i i 
determine 
M = eteeeee - Cc. 
C4 ij tj 


Yn our case M= 4, for cells (1,2) and 1,4), encircled im our cost 
matrix, step 1. We wish now to assign to one cell corresponding tx M, 
as large as possible consistent with the requirement 


that all «,, be non-negative, and that >. x. = 4, ahd Dw, sb 
j ij i aes i 5 


Tf wi increase the assignment in cell (ios 3.) by C7 we must decrease 


7 
by a corresponding amount the assignment in some other cell of row Los in 
order to meet the latter constraints above, This adjustment will cauge an 
unbalance im a different column; so we must continue this compensating 


procedure until we complete a cycle by decreasing by OF the assignment 


im some celi «i, jd 


a 
. ple we fre two candidates for the assignment ce cells 
(j 2 3 (tf rel rhitreé Aan aot 3% ass ae } R: “. 
t,o? Me yl 4). we arbitrarily choose cell (1,2) and formally enter 1 
therein becarvre . mot yet 7 itS value. Making the required ad- 


justmencs we leive che assignment matrix in the following interim form, 
the cycie we have followed, alternately adding and subtracting O.. Ls 


il Tné Subscript L refers to the first iteration, in subsequent 


iterations we use Oe, O neces until the problem is solved. 
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Im acu iteration, at this point there will always be ar least ture 


epi r lee, [ x a? ©]. We make © as large as possible, reaping £4) 
Lj 


oa 
—w 


a) 
Reo 


fp omwrder Lor 4 -.. we 


i | " 
te 


x 2 i in the present example QO, 


non-negative, Hence we set © 2 10 and get the rew assienment macrix 


shown below, 





Assigmient Matrix (step 2) 
ee 
ie ae —_ : Ta. -—— = 
—— or ee 
| | i | | 
a SG : {Oo | 69: 
} | 
! | 4a} 
| ‘ 
ae +9 
| 120 RS OMe 
ee ae FT 
: | ' 
ye | | 
/ /40\ fo 30742 | 
J Aon => 1 me b. : 
TMs Saciuie 7 recuited in a new, basic feasible selution We have 
mae +. PAwSitive and Xi, 7 0. If, at any stage, the introduttion 
of i-})o%.7° £o one route causes two or more positive He to vanish, 
’ : é | 


oJ 


GQ0ELiv@ey +2 undicated. The inclusion of the € perturbation, previous- 


by described, will eliminate the degeneracy. 
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elow wh xe@write Our |Ssipnment matrix (step ¢£) and it@ aseociated 


Stec 2 


Assizgumiat Matta 





we fave @pain encircled two possible routes for the introduction of OF 
Again oir choice is somewhat arbitrary, but route (4,2) has the lessir 
direct cW-t; hence, we choose it. We have indicated the plus and minus 
je ."s above, From the assignment matrix we see ey € 20, in order 


ha 


that X19 not be non-negative. We let OC, = 20, Continuing this 
ra 
process. we go through the following iterations. 


Step 3 
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Step 6 











represento a aogier. 


com 1LOELS =e 
+ 390.1) 


30 


(is 
rei 
a. 


mam (t,., - ©, ph = Di 
ij ij 


ihe cost 1s 


129 + 10(5) + 30.2, + 4013, 


+ 60 + 120 + 40 + 80 + 3H 


2 The Ferd - Polkersen Algorithm. 


Twa matrace¢: are carried along: 


(24 


& 


with "wereore” attached ta the rows and columns. 


(2) An additional matrix of the same size ae C 


Matrix. This matrix is used in making the ascizgnments 


> 
é 


n 


the positive X45 


(1) The cost-weight matrix, This 16 the origina! 


the assilenment matrix 


GUFT * 2ffias 
cose matKeiy C 


aided ttm zere 


B.S. CHT Si uy 


With each row we associate a “surplus” and with 


each column a “shortage; initially these are tne a, 


problem. Censider 


following vableau., 


tie 


iL 


problem previously given which 


we Ghewia the initial cost-weight matrix ard the 


nmer. Preparé tw instrices 


yac matrix C= Pa 30 Leave 
3 


ary i < of tre ori¢ 


‘ 


ab, 


a 
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summarized im ¢ 
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s'4n O0 rew of the cost matrix, pick Bet te® leee enccrvegre SarBq 
ifs Hazee rvs the weieht for that row,  q(Yhis Sten fias been Oa Lana 
Next plece circl#s in the zero matrix correspending to gach appeersrcs of 
these igast values in the cost matrix. (This stew has kis) Leen doe 
ature). 

Then wich each column in the zero matrix that contac! 2 ire? 
ASScciate 4 weight of zero and enter these weights in appripriale posi- 
Custic iz a tow above the cost matrix, Fer each tolumt pur qassastiun: & 
4 iz ‘n three in the example, see matriy beliw) tonoure each entry 
AME 2% fut. Yake the minimum of tnBe®@ reselee ana @nter its 
tis eolumn weight for the column. Ester corrlos am the mer 
MALLE RX & ti fells corresponding to thes. manimal wal Wit? & @hiny 
attecr@d 2. wach row and each column) our Bost mat rir ; =: cee- 
fesbed af @ erSt-wetehe matrix, It wili hersafter te &) cesugmaced xz- 
eiits cf § reecung sien are shown below firr cur example. 


ur 





og ws 


Lo 30 





eo 7 





ih? next scep may be thoughe of as partially satisfying shortages 
in che columns from surpluses in the rows, using only circled positions 
a® mogeinle shipping routes. This step is accomplished as follows: 
seum ch rirst row for the first circle appearing in it. Letermine the 
M.Qimun of the surplus for that row and the shortage for that column, 


Enter this number in the circle and simultaneously reduce the shortage 


arma “urplu? by that amount. Then continue to the second circle aad se 


? 


Cn Jpan the first row is completed go to the seconde. Repear until ali 


e 


rows ave been scanned. Applying this procedure to the zero matrix of our 
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(ev Ve Serereece en to@racin prqevedyrs Which w.1l ewencueley pra- 
vtec es Le  pYaQcedune ebarts wit! $40@i*s) Ura ye %.  - 
mse 10 Ee ceskit rough or piu Eig gue he Peceniae ton 
eat xs we wl VL bt discussed b@icw oct far f@i -s gebotiod the 
bebe lige prow 
Star; Oy iSbeling the rows which stiil neve curpl umes (Yo War exemple 
there LF tciv owe ves thé label (2 foe seach Wem rou, where 0 
15 “2 on f xrplus in that raw, o. refers ta the column of syr- 
p *~ Rm example thiS sifep ghens the following diagram 
fhe 
! to / 
er@aeitar che re 2s ¢€ wilt be ciferrea tc as the ‘‘pSiential Tiow” 
a$8Cc late Li ii 
MY xt we igeg Mmitially, cime cf them are Llsbeled arc tris 
point. We tame @¢acn 'Moeied row and tweep along it booking for circles 
in currentl: wnsabeled cuiumns. Lf such a carche 18 found Label che 
eclumn contaiming it «f ey where £ is the potential flow asseciated 
With the row being anned and R identifies the row. Results of this 


Step ere shown below, 
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wereafter be called the “peten: 


- 


The number £ witil 
iated with the column. 

When the rows have all been scanned, turn to the newly labeied columns. 
Scan each labeled column for circles containing positive entries bet not 
appearing in a previously labeled row. ‘ Any rew in which such a circle lies 
is given a label (2 , or where oF is the column being scanned and 
HX is the smaller of the potential flow of the column and the entry in the 
eircle. Column scanring differs from row scanning in two ways. The circles 
nmust have positive entries, and & is computed differently. Applying this 


sten to our example, yields the following diagram: 








Having scanned all newly labeled columns we turn again to the rows 
scanning the newly labeled ones as before. Continue this procedure 
scanning newly labeled rows and columns until cirher 9 oreakthrough 15 
achieved or no new labeling is possible. When scanning row #3 we achiev- 


ed a breakthrough (i e we labeled a column that had a shortage}. At chat 


point the matzix was labeled as bel w. 





(40 he \30 Qe’ alta aR) 


Breakthrough 


We stop labeling when a breakthrough ts achieved and make adjust- 
ments in the allocations. Suppose we have labeled a column which has a 
shortage, s, with the label (Cf, R.)- In our case s = 60 and the iabel 
is (30, R,). This means that we can increase the flow to the Labeled 
column by an amount h = min (s,£] o The following procedure - used, 

(1) Decrease the shortage by an mount h. (39 in our case) 


(2) Look at the label (£,R. ? avd imcrease the entry in the circle 
where row i and the labeled column intersect by h units. 


The label on row i designates 4 certain column, j. Decrease 
the entry in che circle where row i meets column j by h units. 


a, 
Lo 
Ree? 
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(4) Obtain a new row from the label on column j and increase) 
continue alternately increasing and decreasing the amounts in 
circles by h umits until a row is reacred whose label 14 

7 ; 
tet, bc). 
. Lo af 

(5) When the indicated surplus is reduced by th units, this process 

Ls completed. 


Applied to our example these steps produce the following matrix. 


Arrows indicate the "path" we followed in the adjustment process. 





Tf, at this point all shortages and surpluses are zero, the 
problem is solved. Otherwise, remove all previously applied labels and 
commence the labeling process anew. 

Next time around, we exhaust all possible labelings without achiev- 


ing further breakthroughs and have the matrix shown below. 
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' 10, Ry 39 Ra (SO,Ry “| Gi) 


Although we have labeled all excepre one row and one colunm, we 
have not achieved breakthrough, At this point we turn sur artceantlon 
to the cost-weight matrix which provides 4 means of proceeding. 
Non-breakthrough 

We now change the weights on the cott-welent matrix by sultracting 
k units from the weights of the labeled rows and adding k units to the 
weights of the labeled columns. The value of k is determined by taking 
it as large as possible subjece to the comstraimt c,, + Ww. + w, et) 


= J 


where c. is the ij entry in the cost matrix and w 


‘ is the row weight, 
¢ is the column weight. The determination of k can be accomplished 
as follows 

Take up the cost-weight matrix but igmore for the time being all 


labeled columns shading, covering, or some other distinguishing mark- 


entry in each label- 
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ing is recommended for the labeled columns. Fs 
ed row, compute the algebraic sum of it, its row weight, and its column 


weight. The minimum value obtained this way is k. Enter circles in the 
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r gage All wolwmns have been labeied excep: seed con- 


Sider Owiy pO*itions «1,3%y .2,3), and {» 3). New crscliet are placed in 
pesiticns 2 3! and ¢€4.3), indicated temperarily by diaronds in the zero 
matrix. Subtract & from each labeled row's weight. and add k to eacn 


wan’s weight. Now looting at the zerc matrix; we ipnore all 
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Jabelec co 


lahelea rews and remove all circles réemareing in tabeiem columns oaly 
toe circle in po om (3,2) £8 to cemov@d im cur case. Te has been crmss- 
Si out Glove. 2 check on the procédur€. this proces#® snould never 
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lead th che r sal of a circles with 4 po@itive wniry 


ove procedures provide uS a4 new zerg matrix with which we re- 
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commence the labeling process. A new cost-weLrent satrix is alse shown 
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With very little labeling we chieve breaethrough into column 3. 





This breakthrough provides us the needed 30 omits flow to column 3. and 


the problem is solved. All shortages and surpluses aré now zero. 
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This }5 the same cust as obtained by the 


The above description is reully more unvolved than the actual ase c 
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© of CnSS mentmay sSbitwer is elem as fol Wwews 
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bye + EoD; 43 + @ ee cd] 


other methods a8 act must be 


tt, 


the algorithm Once the ground rules @re underscocd, the manipulation of 
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the matrices is not difficult 
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The authors of this algorithm strate tuey 


compleced a 20 x 20 cptimal assignment problem im about 30 minutes of 


hand computation. The same problem by the simptew method, reqeired weli 


Over an hour. 


3. The Munkres Algorithm, 


Because of its basic similarity to the Ford - Fulkerson method. 


this algorithm wili be stated completely. and then the example will p< 


worked out. Munkres statement of the problem is given to speficy his 


notation. 


Let ; 


D [a J be an mx mn matrix of mon-negative integers, 
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vr. (7 = 1,2,...m) and 
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integers such that Zo 
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=— 


1,2,...m) be positive 


= Cola. 
» J 


Determine values of *1i which minimize the sum Ped, subject to 
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The Exsilowioz statetent Of the A@lgoritior i® Cave. directly from 
Mitered pag: "Siger ltitwe for the Assignee @ a JPewnept «cet tie Pls fete 

We werk with the cost - 7 - [4 7 Tr tte cole se wet the 

probiem we wiil avec Lipa sees aew af the matriy callin, 


if 
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them aovered tines, and we will dist it c8rtaln zerc slemerta 
of the matrix by means of asterisks ore bo oe Tn addition we 
Will a@slgen to e&en elemeprt af the merle 4 atinevesketive auntea x 
which w&) be chief Fp ci AX ~~ PPL Fact elemiar’J 
of Che Natiinx wnhoee sunhs Bsr Teves weit &: Siked 35 mt .o0 Caer 
el.memts wild always 52 verte Bi ditry ©: We Sf the prvb@err wi he 
number c, -2 x will be calied tae (\§ereweesty of the jreh 
| | jt om ee <S ae 
cc Lem atJtnat Staze, and tre wumbes + = ‘ is the discrepancy 
of tne ith rew. These discrepanci 9 (i af @dy $7 3e non-negative ; 
when they ali vanish, the correspredimy Ghote’ are a sslurticr to the 
problem, 
We shalt construct a diagenal in sacnm eral VI i: “YLx and place 
the od, below and the corresponding 2 . abow! toe agondi. Inirctal 
3 tJ 
discrepancies (the ry and c,) will be written to the lefc and above the 
: 3 
matrix. Discrepancies at succeeding etaren will oO sma within dorted 
lines with row discrepancies to th® rvieo® bw -25min dpecrepancies below 
the cost matrix We will demote row ¢1s_cepancties 9. A ,y columa dis- 
crepancies, as fA, i = 1, m j= Epo eee 7 Ld inmmecessary 
I 
clutter zero quotes will not be showr explicitly, bor will be implied by 
a blank above the diagonal im each cell where 4 2er. wusta is assigned. 
Preliminaries 
All quotas are zero, no limes are covered, no geres are Starred or 


primed, 


subtract from each row of the matrix D aut<s smallest emtry, then 


subtract from each column of the resulting matrix its emaklest entry. 


Find a zero Z in the matrix. Tf the discrepancies of both its 


row and its column are positive, increase the guota assigned to Z 


yf 


Ja 


Ke! 
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until the smaller the two discrepanci reduced to zera Repeat 





for each zero in the m&crix. Then cover every ©i3)l@mr wroee discrepancy 


Step 1, Choose a non-covered zero and prime rt. Mwsaiagr tte row tontain- 
ing ic. If che discrepancy cf this row 1s p28iatuve go at once to step 2. 
Otherwise (1f the row discrepancies is zero) comer tha» row, Then scar each 


column. Rée- 


aj 


cwice coveréd essential zero Z2 in the row acd _ncover Zz 
peat umtii1 all zeros are covered, Go to step 3. 


Step 2. tnere if 4a sequence of alternately starred and primed zeros con- 


Structed as follows: Let Z. dencte che umcvvera? O° (there is only one). 
- 
Let ¢, denote the O* in zs column «if any), tet 7, denote the 9° 
a ( 
in Z's rew, Lat 2, cencte the 9* in Z's ‘glum (if any. oimilarly 


continue wnmtil the secuance stops at a 0". Zaye y which has no O* in its 


mo column contains more thaz one OF ana no row more rhan 
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ome 0°, this sequence is unique, It may contain only one 
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Z "Ss row 18 positive, that of Z 
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discrepancies & cols 1s positive, 


and the quota assigned each O* of the sequence Zisoee Lay is positive. 

Let nh be the smallest cf these positive menbers. Increase the quota of 
each O' im the sequence by n, and decrease the jy-ota of each OF in the 
Séquence, by bh. frase all asterisks and primes. Uncryer ali rows and 
cover every colum: whose discrepancy is zers. ktturn te st¢o ?}. 

Step 3. Let k dénsre the smallest non-covered clement of the matrix; it 
wili be positive. Add « ta avery covered row and subtract k from every 
uncovered coiumn. Return to Step 1, without altering any asterisks, primes, 
or covered lines, 


Ovr sample problem will be stated in Munkres nctation and eclyed 


according to his algorithm 





‘G 





4. 


PrSeLiminaer.¢? 


ee 


at 


Pesults 





63 





of Step 1 


Ve 





a Te 


§ 


YY eae oo wr mw eer 


-—-— —- 2 we ow we fF 


30 


Pr 


' 
i 
t 


. 


rr 
Ny 
©: 
Or 


a8) 


th 


rm 
4 


re) a ¥5 


< 


Cur 


Ze Yi 


we 


7 
iy 


ae 





CALcr .€ 
hecame 
2 with 


of Bren 3 


Our 


’ 


« 


atcrix only in céade 


We prime 


n* 


- 


oo” 


a 


Zexre 





The cost 18 G30. as beétor: 


The Similarities between the 


*he manner in which 


ta. oA wy * 
ft K ae ree BS aes 


Js 


Lit 


latter two 


ns dewelcped 


aleoxrithms 1S apparent in 
By ] 





BILLiOyharry 
Dart airenc. EF “Applicaticn of the Sumplex Methsa to a Transport a 
tion Problem", Activity Analysis of Productien and Allocation, Jotm 
7 au —Syee I ES A ee a eg ee ee ee ae ee oe ee 
Wiley ws sons, Tnct., New York, 1951 








C&S, Saul i., Linear Programming Methods and Applications, McCraw- 
Hill Boce Co., Inc., New York, 1958, 





Ip 


APPENT:IX §& 
SOLUTIONS OF A MAXIMAL FiiGW PP IRiIM 
1 ceneral 


In tnis appendix we shall present a smali pacicitated merwork 


{> 


and derermine the maxinial flow thereim by using che tecmisgques described 
ie Chapuer ITi fe s.eil work the problem frrat by che flsoding technigue 
and second by using the minimal-cut procedure, fn each instance, we shall 
restate the essential features of the technique enpioyed, 


2, The Problem. 


—= 


25 


20 


\ . 





For the network WN shawn above determine the maximum tLeady-sState flow, 


From origin 4 to terminal b. 


5B TRO WY Aver. Piooding Te®hn: 
StharShime Ve the origin, assign SurSitdim. Csbbbkic flaws to ali arcs 
Seaving the origin fo saturate thar Vid 1.) 88 Wee © f jueection points as 


new origins and, startup wath the ome cobyect to the greatest capacity 


constraint, scheduié units whenever possible in the following order: 


> 


a 6 N Forwa rd a J jJunee Lon pe VEN Be 2 eee a 3% ni od pase OUC RO ing Arcs : 
‘a fie = . ” . a is *« Pe ; ame 
‘Lateraily - 12 other poinrs of the fer, 
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3 'Batctlenecked’ - if units are left over after steps 1 and #, 
hoe. if ALL Sun@oine and lateral arcs @/™ saturated 


Continue the above d@rocedure Until thé EcGe ible l tiow Cavete vee 
complete retwork and reaches the tecninal 

Eliminate bottlenecks by returning ail exce%* units to the origin 
The validity of the solution can be checker sy inspectidn. If a maximal 
flow has been achieved, there will be ne conmtimuous unsaturated path ex- 


tending from the grigin to the terminal. 





4s Ws urior by Polydyreff s Flooding Techniqye 
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The above represents the res lt® of ile first two srages of tivoding 
Bottlenecks: 6 at function ¢; iv 9t junction d, 
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Ferthear application yrelds: the didaprar bel (Peteh= clirupl vresesrares 
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LO 
One Daieictonat unit 13 bottlenecked at point £F. 

Return bottlenecked units to origin, and the aaximeum flow is seen to 
be 28 units. 
3. The “ithinui-Cut Frocedure, 

Let T be tne chain joining a4 and b whicit 15 tupmoast pi bh Impo 52 
as laree « flow as pessible (T, ky en this chain, thereby saturating 
one of it. arcs. 

Subtract Ke from the capacity of each are of T. 

Remove the previously saturated arc whose capacity 1s now zere Record 
ky Continue this procedure, Eventually, the graph will disconnect and 


the maximal flow 15 established as: 


Sng “ \g 1 
= | 
F = @ t, Lae 
be 
rahe 


where > ke is the amount of flow necessary to saturate the topmost 
ae 


war. th ; 
chain in the i step, and mn is the number of steps required 


for the graph to disconect. 
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6 Solurion by the Minimum-Cot Procedure 


Scep 1 


qT, = <a,c,eée,h,b,) kK, 
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T a rg ane es a ; 


Arcs ¢j and hb caén no lenoer contr bs te to any chain flow from 4 to b 








Step 6 


T. = (a,d.f,g,j3,>) 


Step 3 





Te = (a,d,a@,j,0) 











step $ 
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Step 10 
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On completion of this stap the graph disccnne ts, 


The capacity of N is given As. 


Cap (N) = - 


x 26 wnits. 


ZK = St 1+3+2+ 34242424342 
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APPENDIX C 
GRAPH THEORY 
In Chapters II, TII, and IV of this thesis we have used the lérguage 
ef the mathematical theory of graphs. We include here detinitions and 
theorems from this theory, in order to amplify and clarify some of sur 
statements isi the text. 
1. DEFINITIONS. 


A graph 7 ts a finite, one-dimensional complex, consisting of 
b,c,-....@ and arcs GA (ab), Plac)eccs, d (ce). 


An arc CQ (ab) joins it end vertices a and b; it passes through 
no other vertices and intersects no other arcs ¢xcept in verkice:. 


A chair (path) is 4 set of distinct ares of G which can be arrfnged 
as@ ab), 6 (bc),..., dS (gh) where the vertices a,b,c,... h ®ce distinct, 
A chain does not intersect itself and it joins 1t$s end vertices, 4 and h. 


A cycle is a set of distinct arcs that can be ordered as QA‘ab}, 
Za ey an Oo (ef), oS (fa), the vertices being distinct 4s in the case 
of a chain. 


A graph is connected if each pair of vertices is joined by a chain 





A forest is a graph containing no cycles, and a tree is 4 connected 
forest . 


We may associate with each arc of a graph a positive number called 
its capacity. 


eraph G together with the capacities of its undividua: arcs, is 
- ING network 


a te = ee eee 


origgted or directed graph is a graph in which each arc has a 


1 , ! . 
dixé=pion soerified. In such a graph. CA (aby; and OQ (ba}, if they both 
exr9r mc wabtincte. 





7 Al ¢ han} 


ii some applications with oriented graphs it is desirable to dis- 

tinguich cértain subsets of vertices as origins and certain others as 

destin®rions Other names used in this connection are origins and ter- 
rmals and sources and sinks. The word terminals is also often used 

t> J'stinguish a subset ¢fF nodes in a network, the remaining nodes being 

jen:fion poincs. 





A two-terminal rectwork is a metwork in whith tw princes are rented 
and distinguished botr from eaca other ang ivom the rEdepeine poinr:. 
Usually the first cf these points i¢ called the saute e & 
called the sink. 





A chain flow from 4 source 4 to sink 6 15 4 Crupia ‘¢ e) "ommEbed 
of a chain joining a and b and a non-negative number « represen sige ihe 
flow along C froma to b. 


A network-flow is a collection of chain flows whith has the peoperi y 
that the sum of the numbers k of all chain flows that contain any are is 
mo greater than the capacity of thar are. Ir equality nheides the are is 
Said to be saturated by the flow, 


The flow in a two-terminal network is said to be uni-directional if 
among all the chains from source a to sink b ail chain flows occur 1m the 


same direction in any one arc. 


2 DISCUSSION 


It is often desirable to associate Aan abstract graph with a fo. 
logical graph fi] . This can be doné by representing each vernaxy ji nece. 


as a point in Euclidean three space and @¢ach are 48 a curve jotning tthe 

end vertices. In keeping with the definitions above, these curyet de pot 

intersect except at vertices. A graph is said to be plarar .f the +e cist- 
: 3 am a ; | P 

ed topological graph is planar. Whitney fil gives other criteria f3- 

determining whether a graph is planar, These criteria are stated 10 the 

following two theorems. 

Theorem 1. A graph is planar if, and only if, it has a dval 

Theorem 2, A graph is planar if, and only if, 1t comtaina neit*er +f tis 

following two graphs as subgraphs: 

(A) G' consisting of five nodes a,b,c,d e, and arcs conné 
of these 10des to every other one by a singie are or a Suspended 
chain, (Note. A suspended chain ts a chain in which omly rt 
vertices are met by more than twe arcs. } 

(B) G°° consisting cf two sets of three nodes each a.t,c, and dd.) t 


and arcs connecting each of the nodes of the first set to ewery rode 
of the second set by a4 Sirglée are or a Suspended chain. 


Ll, A topological graph is planar if it cam be mapped b-i “nto a olan: 2? 
a sphere. 
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An obvious représertatic: 
are used to represent nodes and lines to represent ar: 


be shown as numbers adjacent to the arcs t2 whic» the 


tlon 1s usudily denoted by acrows, 


Another represenration cf a greph 18 4 quer matils 


1s a Zraph 18 & CLAGMSse 


the rows and the columns represent nodes. If wa desire 


graph G by a matrix L = Xk , then tye i li 


ij i 


= 


im the grapn and .2 zero ctherwise, We ctuld mocity chig mm 


capacities by replecing the I"s with numbers wileh tepresen 


In an oriented graph iy may differ from oar - Br t 


positive and the other zere, 


aS 


(Le 


Another matrix representation ef 4a graph 145 an ws iden: 


this matrix the nodes are represented by riowe, the arte. 


say an arc is incident on a node (and vice-wersa, if tF 


the node. In an oriented (directed) graph we say thal 


3 


in to 4 node is positively incident on that mode, an 44. 


from a node is negatively imcident on it Tf we repr 
cidence by a +1, negative incidence by a -1, ape ne 
then in the incidence matrix of a direceed grapk eMc* 


+1 and one -1. Rows which represent origins would 
row of exclusively +1"5 would tndicare a carmimal 


be represented by rows containing l’s Gf mixea signe. 
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